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Abstract 

We consider the three-loop corrections to the static potential which are induced 
by a closed fermion loop. For the reduction of the occurring integrals a combination 
of the Grobner and Laporta algorithm has been used and the evaluation of the 
master integrals has been performed with the help of the Mellin-Barnes technique. 
The fermionic three-loop corrections amount to 2% of the tree-level result for top 
quarks, 8% for bottom quarks and 27% for the charm quark system. 

PACS numbers: 12.38.Bx, 14.65.Dw, 14.65.Fy, 14.65.Ha 



1 Introduction 



The potential between a heavy quark and its anti-quark is a crucial quantity both for 
understanding fundamental properties of QCD, such as confinement, and for describing 
the rich phenomenology of heavy quarkonia [1] (see also Ref. [2] for a recent review about 
the static potential). 

Within perturbation theory the static potential can be computed as an expansion in the 
strong coupling as and the inverse heavy-quark mass or, equivalently, in the heavy-quark 
velocity v. The leading order result in v is known up to the two-loop approximation [3-8] 
which has been completed about ten years ago. For the three-loop corrections there are 
only estimates relying on Fade approximations [9] or based on renormalon studies [10]. 

A new feature of the three-loop corrections is the appearance of an infrared divergence 
which was discussed for the first time in Ref. [11]. A quantitative analysis of this effect 
can be found in Ref. [12] (see also Ref. [13]) where a proper definition of the static 
potential within perturbation theory is provided. Furthermore, it is argued that the 
infrared singularities cancel in physical quantities after including the contribution where 
so-called ultra-soft gluons interact with the heavy quark anti-quark bound state (see also, 
e.g., Refs. [14-16]). Higher order logarithmic contributions to the infrared behaviour of 
the static potential have been considered in Refs. [17, 18]. 

In this paper we compute the fermionic contribution to the three-loop static potential 
which is infrared safe. Partial results have already been published in Refs. [19,20]. 

The static potential enters as a building block in a variety of physical quantities. Often at 
three-loop order only estimations are used or the three-loop coefficient — usually called 
as — appears as a parameter in the final result. Let us in this context mention the 
determination of the bottom and top quark mass from the ground state energy of the 
heavy quark system which has been computed to third order in Ref. [16]. The error on 
the mass values due to the unknown three- loop coefficient amounts to 14% (13%) of the 
total uncertainty for the bottom (top) quark. Similarly, 03 enters the calculation of the 
total cross section for top quark threshold production at next-to-next-to-next-to leading 
order (see, e.g., Ref. [21]). The static energy between two heavy quarks has often been 
used in order to compare pcrturbativc calculations with simulations on the lattice (see, 
e.g., Refs. [22-25]). Also in this context the knowledge of 03 is crucial and is expected 
to lead to a better agreement between the two approaches [24]. Last not least let us 
mention the extraction of the strong coupling from lattice simulations where again the 
static potential plays a crucial role [26, 27] and the knowledge of 03 would be highly 
desirable. 

Let us for completeness mention that the one-loop mass-suppressed corrections to the 
static potential have been evaluated in Refs. [28-33], the two-loop terms in Ref. [8]. 
Light quark mass effects have been considered in Ref. [34]. A collection of all relevant 
contributions needed up to next-to-next-to-next-to- leading order can be found in Ref. [15]. 
Furthermore, the two-loop corrections for the case where the quark and anti-quark form 
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Figure 1: Sample diagrams contributing to the static potential at tree- level, one-, two- and 
three-loop order. In this paper only the fermionic corrections are considered at three-loop 
order which excludes diagrams of type (h). 



an octet state have been evaluated in Ref. [35]. 

The remainder of the paper is organized as follows: In the next Section we present details 
of our calculation. In particular we discuss the various types of Feynman diagrams which 
occur at three-loop order and their contributions to the individual colour factors. In 
Section [3] our results are presented and Section H] contains our conclusions. 



2 Calculation 

In the practical calculation of quantum corrections to the static potential one has to 
consider a heavy quark and its anti-quark which interact via the exchange of gluons. In 
Fig. [1] some sample diagrams up to three-loop order are shown. 

In momentum space the static potential can be cast into the form 
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where as denotes the strong coupling in the MS scheme and explicit results for ai [3,4] 
and a2 [5-8] are given below in Eq. (jlj). The infrared logarithm at order follows 
the conventions of Ref. [15] and the renormalization group logarithms \n{fi'^/q'^) can be 
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recovered with the help of 



vr \ TT / 



+ (^) ^ ^f3'oL' + + + ■ ■ ■ , (2) 

where L = ln(yU^/g^) and the coefficients of the /3 function (see, e.g., Ref. [36]) read 

Po = -( — Ca TpTii 

HO 4 \^ 3 ^ 3 ^ ' 

/^i = ^(^Cl-'^CATFrii-iCFTpni 
16 \ 3 3 

+ Y^F^n^ . (3) 

Here, Ca = Nc and Cp = {N^ — 1)1 (2Nc) are the eigenvalues of the quadratic Casimir 
operators of the adjoint and fundamental representations of the SU{Nc) colour gauge 
group, respectively, Tp = 1/2 is the index of the fundamental representation, and rii is 
the number of light-quark flavours. Let us at this point only mention that Oi has the 
colour structures Ca and Tprii where the latter contribution originates from one-loop 
fermionic corrections to the gluon propagator. Note that there is no colour factor Cp 
since this contribution is generated via an iteration of the tree-level result. Consider, e.g., 
the one-loop planar-ladder and the crossed-ladder diagram of Fig. [T](b) with the colour 
factors Cp and Cp — CaCp/2, respectively. It is easy to see (see, e.g., Ref. [37]) that the 
Cp term can be generated by iterations of the leading order diagram in Fig. [U^a) leaving 
only the CaCp term as genuine one-loop contribution. 

Similarly, at two-loop order there are the colour factors C\, CaTpUi, CpTpUi, and (Tprii)^ 
where the latter two originate from loop corrections to the gluon propagator connecting 
the quark and the anti-quark. a2 contains no Cp and CpCa terms since their contribution 
is again generated by iterations of lower-order results. 

The rule that a colour factor Cp can only arise from corrections to a fermion bubble in 
a gluon line also holds at three-loop level. This requires a careful analysis of the colour 
factors for each class of diagram. For example, the colour factor of the graph in Fig. [T](e) 
receives a contributioE0 (C^?— C^/2)Tpn; from two-loop fermionic subdiagram and a factor 
{Cp — Ca/'2') from the remaining crossed box structure. Whereas the complete first factor 
has to be taken into account only the Ca term of the second factor contributes to the 
potential. In a similar manner all diagrams have to be analyzed which leads to the colour 



^In addition to the factor Cp already present in Eq. ([T]). 
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Figure 2: One-, two- and three-loop diagrams. The sohd hne stands for massless rela- 
tivistic propagators and the zigzag hne represents static propagators. 

structures C^, C^TpUi, CaCfTpUi, CpTpUi, C AiTpriiY , CpiTpriiY and (Tprii)^. In this 
paper we compute all coefficients except the one of C^. The results for the structures 
CA{Tpniy, Cp{Tpnif and {Tpriif can be found in Refs. [19,20]. 

At three-loop order there is a new class of diagrams containing a "light-by-light" subdi- 
agram (see Fig. [l](f) for a sample graph) which develops the colour factors d'^'^'^d'^'^'^/NA 
and C\Tpni. Note that these contributions are not connected to iterations and are thus 
already present in QED (i.e. for Ca = 0, df'''^ = I, Na = I and Tp = l). 

Since we use non-relativistic QCD as a starting point for the evaluation of the Feynman 
diagrams the momentum transfer between the quark and the anti-quark represents the 
only relevant scale in the problem. Thus all integrals can be mapped to the two-point 
functions which are shown in Fig. [2] in diagrammatical form. Next to purely massless 
lines originating from the gluon, ghost and light-quark propagators also static lines from 
the heavy quarks are present. The one- and two- loop diagrams have been extensively 
studied in Ref . [37-39] . As far as the three- loop diagrams are concerned one can perform 
a partial fractioning in those cases where three static lines meet at a vertex. This leads to 
many different three-loop graphs involving, however, at most three static lines. Thus any 
resulting integral is labeled by twelve indices one of which corresponds to an irreducible 
numerator. 

Altogether we have to consider about 70 000 integrals (allowing for a general QCD gauge 
parameter ^) which can all be mapped to one of the diagrams shown in Fig. [31 Thereby 
the linear propagators can appear in two variants: either in the form {—v ■ k — iO) or 
(— f ■ k + iO). If the loop momenta, k, I and r, in the upper row of Fig. [3] are chosen as 
the momenta of the three upper lines, then the first diagram appears in two ways: either 
with the product 

{-v-k- iOy^-v ■ I - iOy^^'i-v ■ r - zO)"''", 
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Figure 3: Three-loop diagrams of "ladder" (first three diagrams), "non-planar" (forth 
diagram in upper row) and "mercedes" type (lower row) which have to be considered for 
the fermionic part of a^. 

or with the product 

{-v-k + zO)~"9(-i; ■ / - iQiY"-^'\-v ■ r - iO)"''". 

The second diagram appears with similar propagators where the momenta {/c, Z, r} (in the 
first variant with — iO in all three terms) are replaced by {fc, k — /,r}. The third diagram 
in the upper row of Fig. [3] corresponds to {fc, 1,1 — r} and the fourth one to {fc, /, r}. In 
the case of the "mercedes" type diagrams in the lower row of Fig. [3] one chooses the loop 
momenta fc, / and r as the momenta of the three lower lines. Then the five diagrams 
appear with static propagators of the form {—v ■ k — iO) with momenta {k,k — r,l}, 
{r, k — l,r — I}, {k, r,k — I}, {k, r, /}, {k, r — I, I}, respectively. 

For the calculation of the diagrams we proceed in the following way: They are generated 
with QGRAF [40] and further processed with q2e and exp [41,42] where a mapping to the 
diagrams of Fig. [2] is achieved. In a next step the reduction of the integrals is performed 
with the program package FIRE [43] which implements a combination of the Laporta [44] 
and the Grobner algorithm (see, e.g., Ref. [45]). This leads us to about 100 master 
integrals which have to be evaluated in an expansion in e with the help of the Mellin- 
Barnes technique. Non trivial examples are discussed in Refs. [19, 20] where also explicit 
results are given. We managed to compute all but four coefficients of the e expansion 
analytically. As a crucial tool providing very important numerical cross checks of the 
analytical results we applied the program FIESTA [46] which is a convenient and efficient 
implementation of the sector decomposition algorithm. Finally, let us mention that we 
evaluate the colour factors with the help of the program color [47]. 

In our calculation we allowed for a general gauge parameter ^ in the gluon propagator 
and checked that ^ drops out in the final result. This constitutes a strong check on the 
correctness of our result. 

In order to obtain a finite result one has to renormalize the strong coupling which we 
perform within the MS scheme. The corresponding renormalization constant can, e.g., be 
found in Ref. [36]. 
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3 Results 



Let us in a first step present the results for tfie one- and two-loop coefficients including 
higher orders in e since these terms are needed for the renormalization procedure. We 
obtain 
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where the one- and two-loop results (in the limit e — > 0) can be found in Refs. [3, 4] 
and [5-8], respectively. In Eq. (Hj) is Riemann's zeta function, with the value C(3) = 
1.202057.... 



The three-loop result can be cast in the form 

as ^ 

where the first three coefficients on the right-hand side read 
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where the SU{Nc) colour factors are given by 

N^-l „ 1 
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3 loop 
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0.2228 
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- 0.2655 
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0.25 


0.1172 
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0.02220 


1.701 ■ 
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Table 1: Radiative corrections to the potential where the tree-level result is nor- 

malized to 1 (cf. Eq. (JED). 



In Eq. ([6]) only the coefficient of d'^^'^d'^'^'^ is affected by the limited numerical precision of 
the four coefficients only known numerically which is indicated by the number in round 
brackets. 

We are now in the position to briefly discuss the numerical effect of the new corrections. 
Inserting the results for ai, 02 and in Eq. ([1]) it takes the form 



Vi\q\) 



47rCpa,(|g|) 
^2 



— (2.5833 - 0.2778ni) 

TT 



+ ( — ] (28.5468 - 4.1471n; + 0.0772^^) 

3/(0) \ 
+ (^^) I ^ - 51.4048n, + 2.9061nf - 0.0214nf j + ■ ■ • 

where the ellipses denote higher order terms and /i^ = g ^ has been chosen in order to 
suppress the infrared logarithm. From Eq. ([8]) one observes that both at one- and two- 
loop order the linear ni term is negative and leads to a screening of the (positive) non-n; 
contribution by an amount of about 50% for ni = 5. Also at three-loop order the linear 
ni term is negative and has a sizeable coefficient. Both for 02 and the nf contribution 
is small; the same is true for the nf term of 03. 

In Tab. [T] we present the one-, two- and three-loop results from the square bracket of 
Eq. ([8]) and choose rii according to the charm, bottom and top quark case. In the second 
column we also provide the numerical value of as corresponding to the soft scale where 
fi ^ mqas (and rriq is the heavy quark mass). It is interesting to note that the three- 
loop corrections computed in this paper lead to corrections which are of the same order 
of magnitude as the two-loop corrections, however, with a different sign. In fact, one 
obtains corrections of about —27%, —8% and —2% for charm, bottom and top quarks, 
respectively. Furthermore, let us mention that the unknown constant has to be of the order 
10^ (and positive) in order to significantly reduce the size of the three-loop corrections. 

Let us now compare our explicit calculation with the predictions based on 
Fade approximation. In Refs. [9] and [10] one can find for 03/4^ the results 
{313, 250, 193, 142, 97.5, 60.1, 30.5} and {292, 227, 168, 116, 72, 37, 12}, respectively, where 
the entries in the list correspond to n; = 0, . . . , = 6. A fit to a cubic polynomial in 
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ni leads to 03/4^ ^ 380.9 - 70.42n, + 2.34n;2 + O.OSnf and 03/4^ ^ 362.0 - 72.17ni + 
2.00nf + 0.17nf, respectively. The comparison to Eq. ([8]) shows that the coefficients have 
the correct sign (except the one of nf which is, however, close to zero) and the correct 
order of magnitude. Let us nevertheless mention that the (numerically big) coefficient of 
the linear ni term deviates by about 50%. 

Finally we want to specify our result for V^(|g |) to QED which describes the potential 
of two heavy leptons in the presence of ni massless leptons. Substituting for the colour 
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where a = a{q'^) is the QED coupling in the MS scheme and after the second equality sign 
ni = 1 has been chosen. This corresponds to a bound state of a muon and an anti-muon 
in the presence of a massless electron pair. The coefficients in Eq. ([9]) are significantly 
smaller as in the case of QCD which results in corrections of the order 10^*^ from the 
three-loop term. 

The terms in Eq. ([9]) originate from corrections to the photon propagator plus the "light- 
by-light" -like diagrams as in Fig. [H^f). Thus for n/ = 1 Vqed can be written in the 
form 
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where the photon polarization function is given by 
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with C(5) = 1.036927. . . and L„ = Xnq'^/ml where nie is the electron mass. In Eqs. (11 01) 
and (fTT!) we have used the fine structure constant a. The three-loop relation to a can be 
found in Ref . [48] . 
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4 Conclusion and outlook 



In this letter we report about the calculation of the fermionic corrections to the static 
potential of a quark and an anti-quark. All occurring integrals are reduced to about 100 
master integrals with the help of the program FIRE. The main result can be found in 
Eqs. (E]) where the three-loop coefficients are given for each occurring colour structure. 
The numerical corrections of the new three-loop terms are quite sizeable when applied to 
the system of two charm, bottom or top quarks. However, for a definite conclusion one 
has to wait for the rii independent three-loop coefficient af\ 

The calculation of af^ is currently in progress. We do not expect any conceptual prob- 
lems. However, significantly more Feynman diagrams contribute which leads to many new 
graphs in addition to those shown in Fig. [31 As a consequence also more master integrals 
have to be evaluated. 
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